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We investigate a prototypical agent-based model, the Naming Game, on two-dimensional ran- 
dom geometric networks. The Naming Game [A. Baronchelli et al., J. Stat. Mech.: Theory Exp. 
(2006) P06014.] is a minimal model, employing local communications that captures the emer- 
gence of shared communication schemes (languages) in a population of autonomous semiotic agents. 
Implementing the Naming Games with local broadcasts on random geometric graphs, serves as 
a model for agreement dynamics in large-scale, autonomously operating wireless sensor networks. 
Further, it captures essential features of the scaling properties of the agreement process for spatially- 
embedded autonomous agents. Among the relevant observables capturing the temporal properties 
of the agreement process, we investigate the cluster-size distribution and the distribution of the 
agreement times, both exhibiting dynamic scaling. We also present results for the case when a 
small density of long-range communication links are added on top of the random geometric graph, 
resulting in a "small- world" -like network and yielding a significantly reduced time to reach global 
agreement. We construct a finite-size scaling analysis for the agreement times in this case. 

PACS numbers: 89.75.Fb, 05.65.+b 



I. INTRODUCTION 

Reaching agreement without global coordination is of 
fundamental interest in large-scale autonomous multi- 
agent systems. In the context of social systems, the ob- 
jective is to understand and predict the emergence of 
large-scale population-level patterns arising from empiri- 
cally supported local interaction rules between individu- 
als (e.g., humans). Examples for such phenomena driven 
by social dynamics include the emergence and the evolu- 
tion of languages [H-Q or opinion formation 

The creation of shared classification schemes in a 
system of artificial and networked autonomous agents 
can also be relevant from a system-design viewpoint, 
e.g., for sensor networks (13 . [lij ]. Envision a scenario 
where mobile or static sensor nodes are deployed in 
a large spatially-extended region and the environment 
is unknown, possibly hostile, the tasks are unforesee- 
able, and the sensor nodes have no prior classification 
scheme/language to communicate regarding detecting 
and sensing objects. Since subsequent efficient operation 
of the sensor network inherently relies on unique object 
identification, the autonomous development of a common 
"language" for all nodes is crucial at the exploration stage 
after network deployment. 

To this end, in this paper we consider and slightly 



modify a simple set of rules, referred to as Language 
or Naming Games (NG), originally proposed in the con- 
text of semiotic dynamics [lol [l6| . Such problems have 
become of technological interest to study how artificial 
agents or robots can invent common classification or tag- 
ging schemes from scratch without human intervention 
Will- The original model [TEEMl was constructed 
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to account for the emergence of shared vocabularies or 
conventions in a community of interacting agents. More 
recently, a simplified version of the NG was proposed and 
studied on various network topologies by Baronchelli et 
al. [lO-lil, and by Dall'Asta et al. [H HI The ad- 
vantage of studying a minimal model is that one can 
gain a deeper understanding of the spontaneous self- 
organization process of networked autonomous agents in 
the context of reaching global agreement, and can extract 
quantitative scaling properties for systems with a large 
number of agents. 

In the context of artificial agents, there are other pos- 
sible scenarios when the NG algorithm, in addition to 
being interesting in its own merit in studying agreement 
dynamics on various networks, can also be particularly 
useful from a system-design viewpoint. That can be 
the case when one does not intend the outcome of the 
agreement process among many agents to be easily pre- 
dictable. The actual process of electing a "leader" or co- 
ordinator among sensor nodes may actually be such a sce- 
nario. The leader must typically be a trusted node, with 
possible responsibilities ranging from routing coordina- 
tion to key distribution [26| . Standard leader election 
(LE) algorithms [27l - [3~D ] are essentially based on finding 
global extremum (e.g., maximum) through local commu- 



nications [13-HH • Thus, the elections can be stolen by 
placing a node in the network with a sufficiently high ID 
(e.g., the largest number allowed by the number repre- 
sentation scheme of the sensor chips.) Along these lines, 
a possible application of the NG algorithm is autonomous 
key creation or selection for encrypted communication in 
a community of sensor nodes. Instead of having a central- 
ized or hierarchial key ma nag ement system with domain 
and area key distributors [26[ , group of sensor nodes can 
elect a key distributor or a security key for secure com- 
munications between group members. 

This work is an expanded version of our preliminary 
results |32j . In addition to showing more detailed and 
extended simulation results, we also study and analyze 
different aspects of the behavior of the model, in par- 
ticular, the probability distribution of the agreement 
times and the cluster-size distribution in the NG on ran- 
dom geometrical graphs (RGGs). Further, we construct 
and present finite-size scaling for the agreement times in 
Small- World H (SW)-connected RGGs. The remain- 
der of the paper is organized as follows. In Sec. II we 
briefly review recent results on the NG on various regular 
and complex networks. In Sec. Ill we define and present 
results on the NG with local broadcast on RGGs, moti- 
vated by communication protocols in sensor networks. In 
Sec. IV we present and discuss results on the NG on SW- 
connected RGGs. Section V concludes our paper with a 
brief summary and outlook. 



II. BACKGROUND AND PRIOR RESULTS ON 
THE NAMING GAME 

In the simplified version of the NG, agents perform 
pairwise games in order to reach agreement on the name 
to assign to a single object. This version of the NG was 
investigated on fully-connected (FC) (also referred to as 
mean- field or homogeneous mixing) [2(| [2l[ , on regular 
[22T |. on small- world (SW) [H, H3|, and on scale- free net- 
works [24|, HE)]. In the FC network, each agent has a 
chance to meet with all others and compare their cur- 
rent local vocabularies (list of "synonyms" ) before up- 
dating them. On regular networks, agents have only 
a limited and fixed number of neighbors on a one-, 
two-, etc., dimensional grid with whom they can inter- 
act/communicate. The communication in both cases is 
"local", in that pairs of agents are selected to interact 
and to update their vocabularies. The basic algorith- 
mic rules of the NG are as follows [2(| H^|. A pair of 
neighboring nodes (as defined by the underlying com- 
munication topologv), a "speaker" and a "listener", are 
chosen at random [35j | . The speaker will transmit a word 
from her list of synonyms to the listener. If the listener 
has this word, the communication is termed "successful" , 
and both players delete all other words, i.e., collapse their 
list of synonyms to this one word. If the listener does not 
have the word transmitted by the speaker, she adds it to 
her list of synonyms without any deletion. 



Among the above rules, the restriction to a single ob- 
ject [2(J,l2jJ strongly reduces the complexity of the model, 
compared to a more general case where the naming pro- 
cess of multiple objects can be performed simultaneously. 
From a linguistic viewpoint, this rather strong restric- 
tion is equivalent to preventing homonymy, and instead, 
treating all objects independently. This strong assump- 
tion can be more realistic for a system of artificial agents, 
where agents assign random numbers (e.g., chosen from 
2 31 integers) as "words" to new objects. In this case, 
the number of potential words can be far grater than the 
number of objects, and the probability that two players 
invent the same word for different objects (hence giving 
rise to homonymy) is negligible. 

It was found that employing the above local rules (pair- 
wise interactions), after some time, the agents vocabu- 
laries converge to a unique vocabulary shared among all 
agents |20l - [23| . The major differences between the NG 
on FC graphs and on regular low-dimensional grids arise 
in the scaling of the the memory needed to develop the 
common language before convergence occurs, and in the 
scaling of the time t c needed to reach global agreement. 
(The memory need in the present context is the typical 
value of the largest number of words an agent may posses 
throughout the evolution of the game [20, H3-) I n the 
FC network, the convergence process to global agreement 
is fast [t c ~ 0(N 1 / 2 ) for N agents], but large memory 
[©(./V 1 / 2 )] is needed per agent [20|- For a regular two- 
dimensional network (or grid) , spontaneous evolution to- 
ward a shared dictionary is slow [t c ~ O(N)], but the 
memory requirement per agent is much less severe [0(1)] 
22j. When the NG is implemented on Watts-Strogatz 
33| SW networks, the agreement dynamics performs op- 
timally in the sense that the memory needed is small, 
while the convergence process is much faster than on the 
regular networks [t c ~ O(N 0A ), closer to that of the FC 
network] p3| . 

Sensor networks, which are motivating our study, are 
both spatial and random. As a large number of sensor 
nodes are deployed, e.g., from vehicles or aircrafts, they 
are essentially scattered randomly across large spatially- 
extended regions. In the corresponding abstract graph, 
two nodes are connected if they mutually fall within each 
others transmission range, depending on the emitting 
power, the attenuation function and the required min- 
imum signal to noise ratio. Random geometric graphs 
(RGGs), also referred to as spatial Poisson/Boolean 
graphs, capturing the above scenario, are a common and 
well established starting point to study the structural 
properties of sensor network, directly related to coverage, 
connectivity, and interference. Further, most structural 
properties of these networks are discussed in the litera- 
ture in the context of continuum percolation (36M38| . 

The common design challenge of these networks is to 
find the optimal connectivity for the nodes: If the con- 
nectivity of the nodes is too low, the coverage is poor 
and sporadic. If the node connectivity is too high, in- 
terference effects will dominate and result in degraded 
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signal reception [39|-[43l|. From a topological viewpoint, 
these networks are, hence, designed to "live" somewhere 
above the percolation threshold. This can be achieved 
by adjusting the density of sensor nodes and controlling 
the emitting power of the nodes; various power-control 
schemes have been studied along these lines [H, 52, S3| • 
In this paper we consider RGGs in two-dimensions above 
the percolation threshold, as minimal models for the un- 
derlying network communication topology. Further, we 
consider RGGs with an added small density of "ran- 
dom" long-range links. The resulting structure resem- 
bles small- world (SW) networks l33l |44| . also well stud- 
ied in the context of artificial [45|, |46| and social systems 
[i3.l47j.The focus of this work is to study the NG algo- 
rithm on these spatially-embedded random graphs. 



III. NAMING GAMES ON RANDOM 
GEOMETRIC NETWORKS 

A. Random Geometric Graphs 

As mentioned above in the Introduction, first we con- 
sider random geometric graphs in two dimensions (36l - [38| 
as the simplest topological structures capturing the es- 
sential features of ad hoc sensor networks. N nodes are 
uniformly random distributed in an LxL spatial area. 
For simplicity we consider identical radio range R for all 
nodes. Two nodes are connected if they fall within each 
other's range. An important parameter in the resulting 
random geometric graph is the average degree k (defined 
as the average number of neighbors per node), k=2K/N, 
where K is the total number of links and N is the num- 
ber of nodes. In random geometrical networks, there is 
a critical value of the average degree, k c , above which 
the largest connected component of the network becomes 
proportional to the total number of nodes (the emergence 
of the giant component) (36i - [38| . For two-dimensional 
RGGs /c c «4.5 [38]. There is a simple relationship be- 
tween the average degree k, the density of nodes p=N/L 2 , 
and the radio range R of the nodes [3a - [38l ]. k = pnR 2 , 
which can be used to control the connectivity (average 
degree) of the network. 

B. The Naming Game with Local Broadcast 

We consider the Naming Game on random geometrical 
graphs. In the original context of the NG, agents try to 
reach agreement in finding a unique "word" for an object 
observed by them. In one of the above proposed potential 
applications, agents try to generate a shared unique key 
for encrypted communication. For simplicity, we will use 
the term "word" for the latter as well when describing 
the algorithm. 

Motivated by communication protocols employed by 
sensor nodes, we modify the communication rules to 
make them applicable for sensor networks. Instead of 



pairwise communications, nodes will initiate broadcast 
(to all neighbors) in a continuous-time asynchronous 
fashion. In this paper we consider the initial condition 
when the "vocabulary" of each node is empty. At every 
elementary time step, a node is chosen randomly out of 
N nodes (mimicking Poisson asynchrony for large N). 
This node (the "speaker") will broadcast a word from 
her list of "synonyms" ; if her list of synonyms is empty, 
the speaker randomly invents a word; if she already has 
several synonyms, it randomly chooses one. Her neigh- 
bors (the "listeners") compare their vocabularies with 
the word transmitted by the speaker. If a listener has 
this word, she considers the communication a success, 
and she deletes all other words, collapsing her list of syn- 
onyms to this one word. If a listener does not have the 
word transmitted by the speaker, she adds it to her list 
of synonyms without any deletion. If at least one listener 
had the word transmitted, the speaker considers it (at 
least a partial) success, and (somewhat optimistically) 
collapses her list of synonyms to this one word. At every 
step, the "success" rate S is defined as the fraction of lis- 
teners who were successful (i.e., those that had the word 
transmitted by the speaker). From the above it is clear 
that one of the successful listeners, if any, has to report 
the outcome of the "word matching" to the speaker. In 
order to achieve that efficiently, in real sensor-network 
implementations one can employ the "lecture hall" algo- 
rithm [3 [49[ . In this paper time t is given in units of 
one "speaker" -initiated broadcast per node. The main 
difference between the above algorithm and the one in 
Refs. [20l - [23l | is the broadcast (instead of pairwise com- 
munications) and the underlying network (RGG in this 
paper) to capture the essential features of the NG in sen- 
sor networks. 

When starting from empty vocabularies, agents in- 
vent words randomly. After time of 0(1) [on average 
one speaker-initiated broadcast per node], 0(N/(k + 1)) 
different words have been created. Following the early- 
time increase of the number of different words Nd(i), 
through local broadcasts, agents slowly reconcile their 
"differences" , and eventually will all share the same word. 
First, a large number of small spatial clusters sharing the 
same word develop. By virtue of the slow coalescence of 
the interfaces separating the clusters, more and more of 
the small clusters are being eliminated, giving rise to the 
emergence of larger clusters, eventually leading to one 
cluster in which all nodes are sharing the same word. As 
suggested by Baronchelli et al. [22[ , this late-time process 
is analogous to coarsening, a well-known phenomenon 
from the theory of domain and phase ordering in physi- 
cal and chemical systems [5(| ■ Figure [1] shows snapshots 
of vocabularies of the nodes at different times. For later 
times, group of nodes which already share the same word, 
slowly coarsen, until eventually only one domain prevails. 
This behavior is also captured by Fig. [2fb), tracing the 
number of different words as a function of time Nd(t), 
eventually reaching global agreement, Nd = 1. 
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FIG. 1: (Color) Snapshots of the time evolution of the con- 
tents of the agents' word lists during the process of reaching 
global agreement on RGG for N = 1,000 nodes at time (a) 
t = 1; (b)_t = 43; (c) t = 169; (d) t = 291. The average 
degree is fc«12. Initially, the word lists are empty for all 
agents. Time, as through the paper, is measured in units of 
"speaker" -initiated broadcasts per node. Different colors cor- 
respond to different words, with black indicating nodes with 
multiple words. After the early-time increase in the number 
of different words in the systems, small spatial clusters shar- 
ing the same word quickly form, then subsequently "coarsen" 
until eventually only one global cluster prevails. 



C. Basic Scaling Considerations and Analogy with 
Coarsening 

Before turning to the detailed discussion of our simu- 
lation results, we first sketch the framework of coarsen- 
ing theory [5Cj|], applicable to the observed late-time dy- 
namics of the NG on regular d-dimensional lattices [22| ■ 
Coarsening has also been observed in other models rel- 
evant to opinion formation and social dynamics [H, l5l| . 
Unlike other minimalist (typically two-state) models of- 
ten employed to study opinion formation ML such as the 
one studied by Sznajd-Weron & Sznajd |9], the Voter 
model [ll], [H^, or the majority rule model [Ij], in the NG, 
each agent can be in an unlimited number of discrete 
states (corresponding to a chosen word). Further, at any 
instant before reaching global consensus, an agent can 
have different possible words for the object. Because of 
the potentially unlimited number of discrete states the 
agents can assume, the late-stage evolution of the NG 
resembles that of infinite-state (Q=oo) Potts model [53l — 

13. 

While RGG is a random structure, it is embedded in 
two dimensions, and we also attempt to employ elemen- 
tary scalin g a rguments from coarsening theory. Accord- 
ing to Ref. [22[, on regular <i-dimensional lattices, the typ- 
ical size of domains (each with already agreed upon one 
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FIG. 2: (Color online) Time evolution of the relevant ob- 
servables in the Naming Game in the fully-connected (FC), 
two-dimensional regular (with four nearest neighbors), and 
random geometric networks (RGG) for iV=1024, averaged 
over 1,000 independent realizations; (a) the total number of 
words in the system N w (t); (b) the number of different words 
Nd(t); (c) the average success rate S(t). The average degree 
of the underlying RGG is fc^l2. Data for the FC and 2d reg- 
ular networks are reproduced by our simulations, following 
Refs. [2(| 52 1 for comparison. 



word) is governed by a single length scale ~ t 1 with 
7=1/2, analogous to that of domain formation in systems 
with a non-conserved order parameter [50j |. Thus, in d 
dimensions the average domain size C{t) follows 



C(t) ~ £ d (i) ~ t d ~< 



(1) 



and the total number of different words N4 at time t 
scales as the typical number of domains 



N d (t) 



N 



N 
t 1 ^ 



(2) 



Further, the total number of words N w (N w /N being the 
average memory load per agent), at this late coarsening 
stage, can be written as the number of nodes N plus 
the number of nodes with more than one (on average, 
between one and two) words, separating the different do- 
mains. It is of order of typical number of domains times 
the typical length of the interface of one domain, yielding 



N w (t) - N 



N 



N 



N 



(3) 



Similarly, the "failure rate" for word matching, 1— S(t), 
(where S(t) is the success rate) scales as the fraction of 
nodes at the interfaces separating domains with different 
words 



1 - S(t) ~ ~ — 



(4) 
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The main feature of the above power-law decays (up to 
some system-size dependent cut-offs) is that the number 
of different words Na, the total number of words N w , 
and the success rate S(t) only depend on t through the 
characteristic length scale £(t). Further, for the typical 
time t c to reach global agreement or consensus, one has 
e*(*c)~JV, i.e., 

t c ~ N 1 '^ . (5) 

Unless noted otherwise (as in Sec. D.l and 2), our no- 
tation, t c , Nd(t), N w (t), C(t), and S(t) refer to the 
ensemble-averaged values of these relevant observables. 

D. Simulation Results 

Relevant quantities measured in the simulations are 
the total number of words in the system N w (t) (corre- 
sponding to the total memory used by the agents for word 
allocation at time t), the number of different words Nd(t), 
and the average size of domains/clusters C{t). Figure [2] 
displays the time evolution of these three quantities for 
the RGG, compared to the fully connected (FC) and 
to the 2d regular networks. Here, for the comparison, 
we reproduced the corresponding data of Refs. [2(1 [22j. 
The behavior of the NG on RGG is qualitatively very 
similar to that of the NG on 2d regular graphs. After 
time of 0(1), 0(N/(k + 1)) different words have been 
invented [Fig. HJb)]. N w (t) also reaches its maximum 
in time of 0(1) [Fig. Ht a )]- Focusing on the late-time 
behavior of the systems, plotting N w (t)/N—1, Nd(i)/N, 
and C(t) vs t on log-log scales, confirms the power-law 
decays associated with the underlying coarsening dy- 
namics, predicted by Eqs. (J3J, and ([1]), respectively 
[Fig. |3] . From the data for C(t), we obtain 2 7 =0.79±0.01 
[Fig. Etc)], while from the data for N d (t) and N w (t), 
we extract 2 7 =0.74 ± 0.01 and 7=0.36 ± 0.01, respec- 
tively [Figs. [2tb,a)]. Based on our finite-size results, we 
can only conclude that the coarsening exponent is in the 
range 0.35<7<0.40 for the NG on two dimensional RGG. 
Different exponent values extracted from different ob- 
servables for finite systems long hindered the precise de- 
termination of the coarsening exponent in the closely re- 
lated large- Q Potts model [56H581 ]. There, employing ad- 
vanced Monte Carlo renormalization (MCRG) schemes, 
it was shown that the coarsening exponent (within er- 
ror) is 1/2 (58[. However, finite-size effects and very 
strong transients, in part due to "soft domain walls" and 
domain- wall intersections ( "vertices" ) can produce values 
significantly smaller than 1/2 extracted from standard 
MC methods f56l - [58| . such as ours. 

Measuring the time to global agreement, averaged over 
1,000 independent runs (each on a different RGG net- 
work realization), we also obtained the scaling behavior 
of the agreement time, t c ~N ulo±om and t c ~N lm±om 
for fc«12 and fc«50, respectively, as shown in Fig.@] The 
corresponding scaling exponents both somewhat deviate 
from the one predicted by Eq. ([3]) with the exponent 
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FIG. 3: (Color online) Time evolution of the relevant ob- 
servables in the Naming Game in random geometric networks 
(RGG) for three system sizes on log-log scales, averaged over 
1, 000 independent realizations; The average degree of the un- 
derlying RGGs is fc«12.(a) the normalized total number of 
words in the system N w (t)/N—1; (b) the normalized number 
of different words Nd(t)/N; (c) the average domain size C(t). 
The straight line segments correspond to the best-fit power- 
law decays N w {t)/N-l~t-°- 36 , N d {t)/N~t~°- 74 , C(i)~t ' 79 
for (a), (b), and (c), respectively. 

1/(27). This deviation is possibly due to strong finite-size 
effects, dominating the very late stage of the agreement 
dynamics. 

For RGGs with many nodes, a relevant control pa- 
rameter is the average number of neighbors (or average 
degree) k. For sensor-network-specific implementations, 
as noted earlier, k can be adjusted by increasing either 
the density or the communication range of the nodes. 
We performed simulations of the NG for different aver- 
age neighborhood size k, as shown in Fig. [5] The results 
indicate that the scaling properties (in terms of N) of the 
time evolution of the agreement process do not change. 
The typical convergence times, however, are significantly 
reduced by increasing the neighborhood size. A closer ex- 
amination of the convergence time reveals that, for fixed 
N, it scales as t c ~ fc~ 2 - 6 , in the sparse- network limit 
(k<^N) in two-dimensional RGGs. 

1. Agreement-time distributions 

In addition to the average agreement (or convergence) 
time (i c ) (time until global agreement is reached), we also 
measured the standard deviation At c [Fig. @], and con- 
structed the probability density (normalized histograms) 
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FIG. 4: (Color online) Average and the standard deviation of 
the convergence time t c until global agreement is reached, as 
a function of the number of nodes on log-log scales, averaged 
over 1, 000 independent realizations of the NG on RGG. The 
average degree of the underlying RGGs is fc«12 (squares) and 
k~50 (circles). The straight lines correspond to the best-fit 
power-laws with exponents 1.10, for both the average (solid 
squares) and the standard deviation (open squares) of RGGs 
with fc«12, and 1.07 for those of RGGs with fc«50, respec- 
tively. 
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FIG. 5: (Color online) Time evolution of the (scaled) number 
of different words starting from an "empty" word list initial 
condition, for various average degree k on log-log scales. The 
number of nodes is iV=l,000. The straight line segment in- 
dicates the asymptotic power-law decay as determined earlier 
[Figs. [3jb)], independent of the neighborhood size k. 



2. Cluster-size distributions 



P(t c ,N) for N nodes [Fig. [6]. Since in this subsection, 
we analyze the full probability density of this observable, 
we use brackets for denoting the ensemble-averaged value 
of the convergence time, (i c ), while t c alone denotes the 
stochastic variable, corresponding to a measurement in a 
single realization of the NG. 

Up to the system sizes we could simulate, the standard 
deviation, within error, scales in the same fashion with 
the number of nodes as the average itself, At c ~N 1A0 
(k w 12) and At^N 1 - 07 (k » 50) [Fig.g]. [Suppressing 
large average convergence times and the corresponding 
large standard deviations (through modifying the net- 
work communication topology) will be addressed in the 
next section.] 

Further, the shape of the histograms, for sufficiently 
large systems, remains invariant [Fig. [B]. Thus, intro- 
ducing the scaled convergence time x — t c /(t c (N)), the 
corresponding scaled probability densities p(x) for differ- 
ent system sizes collapse onto the same curve. [Fig. |6]Jb)] . 
The above findings indicates that the convergence-time 
distribution for the NG is governed by a single scale (t c ), 
hence can be written as 

P(tc N) = —L-p( tc /(t c (N))) . (6) 

The distributions exhibit exponential tails for large argu- 
ments [Fig. [Hb) inset] , a characteristic feature of opinion 
dynamics governed by coarsening [66j . 



We also studied the probability distribution of the sizes 
of the clusters during the agreement dynamics P(C, t) 
(the normalized histogram of the sizes of domains with 
different words at a given time) [Fig. El^a)]. Similar to 
the previous subsection, we analyze the full probability 
density of this observable, hence we use brackets for de- 
noting the ensemble-averaged value of the cluster size in 
the system at time t, (C(t)), while C alone denotes the 
stochastic cluster-size variable (sampled at an instant t 
in a single realization of the NG). 

Since the agreement process is governed by coarsen- 
ing, one expects that this distribution exhibits dynamic 
scaling, i.e., 



p (C,t) = -J— p (C/(C(t))) . (7) 



Thus, p{x), the distribution of the scaled cluster sizes 
x = C/(C{t)) remains invariant for different times. 
Our simulations confirm this picture, except for very 
early times (growth phase with initial domains forming) 
and for very late times (where finite-size effects domi- 
nate) [Fig. [7[b)]. The cluster-size distribution exhibit 
exponential- like tails for large arguments [59l - l6ll . W\ , as 
can be seen in Fig. E£b). 
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FIG. 6: (Color online) (a) Probability densities of the conver- 
gence time for three systems sizes. Data are gathered from 
100, 000 independent realizations of the NG on RGG. The av- 
erage degree of the underlying RGGs is fc~50. (b) Probability 
densities for the scaled variable x = t c /{t c (N)) for the same 
data. The inset shows the same scaled histograms on log-lin 
scales. 



FIG. 7: (Color online) (a) Probability densities of the cluster 
size at different time of the agreement dynamics. Data are 
collected through 100, 000 independent realizations of the NG 
on RGG. The system size of the underlying RGGs is TV = 1000 
and the average degree is fc«12. (b) Probability densities for 
the scaled variable x = C/{C(t)} for the same data as in 
(a) on log-lin scales. The solis curve represents the best fit 
exponential-like tail, oc exp(— 1.24a; 1 ' 12 ). 



IV. NAMING GAMES IN 
SMALL- WORLD-CONNECTED RANDOM 
GEOMETRIC NETWORKS 

In light of recent results on NG on one-dimensional SW 
networks [23| . we now consider accelerating the agree- 
ment process by adding random long-range communica- 
tion links between a small fraction of nodes of the RGG. 
Such networks have long been known to speed up the 
spread of local information to global scales [33|, 13, 147L 
|68| , with applications ranging from synchronization prob- 
lems in distributed computing [451 ] to alarm-detection 
schemes in wireless sensor networks (46j . For sensor net- 
works, this can be implemented either by adding a small 
fraction of sensors equipped with long-range unidirec- 
tional antennas ("physical" long-range connections) or 
by establishing designated multi-hop transmission pat- 
terns ("logical" long-range connections) between certain 
nodes (69j . 

We construct the small-world-like RGG (SW RGG) 
as follows. We start with the original RGG (embed- 
ded in d dimensions, where d=2 in this paper). Then 
we add "long-range" links (or "shortcuts" ) between ran- 
domly chosen nodes in such a way that the total number 
of long-range links per node (the density of random links) 
is p. This SW construction differs slightly from the origi- 
nal Watts-Strogatz one [33| (also used by Dall'Asta et al. 
[23|). where random links are introduced by "rewiring" 



some of the original connections. The resulting network, 
however, has the same universal properties in the small- 
p, large- N limit [jfj, which is the center of our interest. 
Further, it is also motivated by actual implementations 
in sensor networks, where long-range "channels" are es- 
tablished in addition to the existing local ones. 



A. Basic Scaling Considerations 

Before presenting simulation results, using scaling ar- 
guments, one can obtain an order of magnitude estimate 
for the crossover time t x present in the SW RGG and 
for the time to reach global agreement t c [H)]. In SW 
networks, embedded in d dimensions, the typical (Eu- 
clidean) distance between nodes with shortcuts scales as 
hw^P -1 ^ [70j - t73 |. Starting from empty initial word 
lists word, for early times (following the creation of 
0(N/(k + l))) different words in the system), the system 
will exhibit coarsening, until the typical linear size of the 
growing domains, £(t)~i 7 , becomes comparable to Isw- 
(Here, both lengthscale measures are understood in terms 
of the underlying Euclidean metric.) After that time, the 
agreement process is governed by the presence of random 
long-range connections, yielding mean-field-like behav- 
ior. Hence the crossover from d-dimensional coarsening 
to mean-field-like dynamics occurs when t*> 1 ^p~ x l d , yield- 
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ing 



P 



(8) 



In a system of N agents, the above crossover is only dis- 
played if the convergence time of the original system with 
no random links would exceed the above crossover time 
N 1 / d ~t^ > p~ 1 / d i ; which is equivalent to the condition for 
the onset of the SW effect N^p- 1 jH, [7(J. Following 
the above system-size independent crossover time, the 
agreement dynamics is of mean-field like, and one can 
expect to observe a scaling behavior closer to that of FC 
networks [2(j. In particular, the time to reach global 
agreement is expected to scale as [23| 



tr 



N l/2 



(9) 



a significant reduction compared to that of the "pure" 
RGG with no long-range links where t c ~iV . 



B. Simulation Results 

Simulating the NG on SW RGGs qualitatively con- 
firms the above scaling scenario. Following the very 
early-time development of 0(Nj(k + 1)) different words, 
the system of SW-networked agents, exhibits slow coars- 
ening, with only small corrections to the behavior of the 
pure RGG [Fig. [8]. In fact, this early-time coarsening 
on SW RGGs is slightly slower compared to pure RGGs 
due to the effective pinning of interfaces near the short- 
cuts |, [Zl [Zi- In the NG on SW networks, how- 
ever, the agreement process only slows down [23], but is 
not halted by "frozen" (metastable) disordered configura- 
tions [l,[zij|. After ap-dependent crossover time [Eq. flSJ], 
(when the typical size of the growing clusters becomes 
comparable to the SW length scale), an exponential con- 
vergence begins to govern the agreement process. This 
final-stage fast approach toward consensus sets in earlier 
for increasing values of the density of shortcuts p, yield- 
ing a significantly reduced convergence time compared to 
that of the NG on the "pure" RGG. The temporal be- 
havior of the relevant observables for various values of p 
can be observed in Fig. [8] 

Plotting the convergence time vs the density of 
long-range links, as shown in Fig. [9ja), suggests that 
(for sufficiently large but fixed N) the convergence 
time approaches an asymptotic power-law t c ^p~ B with 
s=0.79±0.01 [23J. On the other hand, for fixed p and 
increasing N, the convergence time increases with N, 
t c ~N asw , with a sw =0.31±0.01 [Fig.^b)]. The agree- 
ment process is much faster than on a two-dimensional 
regular grid or RGG and is closer to the anticipated 
mean- field-like behavior [Eq. ©] (23J. Thus, in the 
small- world regime (Np^>l) the convergence time de- 
pends on both the system size and density of random 
links, £ c ~iV QSW /p s . 
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FIG. 8: (Color online) Time evolution of the (scaled) (a) to- 
tal number of words, (b) number of different words and (c) 
the average cluster size for SW RGGs on log-log scales, start- 
ing from an "empty" word list initial condition, for various 
density of long-range links p, averaged over 1, 000 indepen- 
dent realizations of the NG on RGG. The number of nodes is 
7V=5, 000 with average degree fc~12. 



1. Finite-size scaling for the agreement time on SW-RGGs 

In the pure- RGG limit (Np<£.l), t c only depends on N, 
t c ~N aRGG with a RGG wl.l0 [Fig.EJb)] (since, essentially 
there are no shortcuts in the system). On the other hand, 
as seen above, in the SW-regime (Np^$>l), the agreement 
time scales as t c ^N asv/ /p s . One then can construct the 
full scaling behavior of t c (p, N), capturing the above two 
finite-size behaviors as limiting cases on SW-connected 
RGGs, 



AfasvJ 

t c (p,N) —f(Np), 

pS 

where f(x) is a scaling function such that 



fix) 



x s if x<l 
const, if x^>l 



(10) 



(11) 



The pure RGG limit (Np-^.1) is recovered, provided that 

t c ~(N asw /p s )(Np) s ^N asv ' +s ^N aRGG , i.e., 



"RGG — «SW + S . 



(12) 



Our measured "phenomenological" exponents 
Q!RGG~1-10, asw~0.3I, and sw0.79, satisfy the above 
proposed asymptotic scaling relation. For analyzing our 
data, Eq. (fTO)) can also be rewritten as 



(Nn) asvj 



1 



-g(Np) , (13) 
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FIG. 9: (Color online) Average convergence time t c for SW 
RGGs. (a) as a function of the density of shortcuts for various 
system sizes. The inset shows the same data on log-log scales. 
The straight lines corresponds to an estimate of the associated 
(asymptotic) power-law. (b) as a function of the number of 
nodes on log-log scales for various density of long-range links 
p. The curves shown are obtained by averaging over 1,000 
independent realizations of of the NG on RGG. The average 
degree of the underlying RGGs is fc~12. 

where g(x) — x QSW f(x). Thus, plotting t c p aRGG vs Np 
should yield data collapse, together with the asymptotic 
small- and large- argument exponents of g (x), oirgg and 
asw, respectively [Fig. ITU]. 

j xasw iix>>1 ■ (14) 



V. SUMMARY AND OUTLOOK 

In this paper, we studied a prototypical agent-based 
model, the Naming Game, on Random Geometric Graphs 
and SW-connected RGGs embedded in two dimensions. 
While the underlying RGG communication topology is 
motivated by large-scale sensor networks, the NG on 
these networks captures fundamental features of agree- 
ment dynamics of spatially-embedded networked agent- 
based systems. We have found that, qualitatively sim- 
ilar to two-dimensional regular networks [22j], the NG 
on RGG can be reasonably well described by the physi- 
cal theory of coarsening. In particular, local clusters of 
nodes sharing the same word quickly form, followed by 
slow coarsening of these clusters in the late stage of the 
dynamics. The typical length scale grows as £(i)~£ 7 with 
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FIG. 10: (Color online) Scaled plot of data shown in Fig. [9] 
as suggested by the finite-size scaling argument [Eq. (|13|>] . 
The straight line segments correspond to the best-fit (asymp- 
totic) power-law behavior of the scaling function g(x) with 
exponents 1.10 and 0.31, for small and large arguments, re- 
spectively, as described in the text [Eq. (|14l) ]. 



the coarsening exponent estimated to be 0.35<7<0.40. 
Our simulation results also indicate that the average 
time to reach global agreement is of 0(j\r 1O8±0 - 03 ) (for 
fixed average degree). The above results imply that, at 
least for the range of finite system sizes studied here (up 
to N — 5,000), the characteristic length scale in two- 
dimensional RGGs grows slower than 1/2. This devia- 
tion, in part, may very well be attributed to the effec- 
tively small system sizes that we could study. Similarly 
strong transients and finite-size corrections, due to the 
presence of "soft domain walls" and "vertices" (domain- 
wall intersections), also made the precise determination 
of the asymptotic coarsening exponent difficult in the 
two-dimensional large-Q (effectively Q=oo) Potts mod- 
els (5314581 ] . On the the other hand, based on our Monte 
Carlo studies, we cannot rule out the possibility that the 
deviation from the 7=1/2 coarsening exponent is the re- 
sult of the inherent local random random structure of 
RGGs (in contrast to regular two-dimensional grids [22j). 

While in this paper we did not address the message 
complexity of the NG explicitly, one can make an order 
of magnitude estimate for the typical number of mes- 
sages needed to reach global agreement on RGGs for an 
efficient implementation. (In sensor networks, this quan- 
tity is also relevant since it corresponds to the global 
energy consumption.) Once the coarsening process be- 
gins, nodes inside the clusters have reached agreement 
with all their neighbors, of which they are readily aware, 
hence, they no longer have to initiate broadcasts any 
longer. Thus, only these "active" nodes, found at the 
interfaces between these cluster (which have at least one 
neighbor with different words) , will initiate broadcast for 
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word matching. Using that the number of nodes at the 
interfaces scales as N/t 1 [Eq. ([3])], and integrating this 
expression up to i c ^iV 108±0 03 , one finds that the total 
number of messages needed to be exchanged until global 
agreement is reached is of Q^N 1 - 68 ^ - 05 ^ 

In an attempt to accelerate the agreement process by 
changing the communication topology between agents, 
we also studied the SW-connected version of the two- 
dimensional RGG. By adding a small density of shortcuts 
"on top" of the RGG, resulting in a SW-like network, 
the convergence time is strongly reduced and becomes of 
O(iV - 31 ), similar to the behavior of NG on the Watts- 
Strogatz SW network [23j , 

In future works we will investigate the NG on more re- 
alistic communication topologies, motivated by and rel- 
evant to wireless sensor networks, in particular, random 
spatial networks with heterogeneous range distribution, 



minimum- node-degree networks [75j . and also networks 
with dynamically changing connectivities. 
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